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Abstract. We report on an XMM-Newton observation of the X-ray afterglow of the Gamma Ray Burst GRB 
011211, originally detected by Beppo-SAX on 11th December 2001. The early afterglow spectrum obtained by 
XMM-Newton, observed 11 hours after the initial burst, appeared to reveal decaying H-like Kct emission lines of 
Mg, Si, S, Ar and Ca, arising in enriched material with an outflow velocity of order 0.1c (Reeves et al. 2002). 
This was attributed to matter ejected from a massive stellar progenitor occurring shortly before the burst itself. 
Here, we present a detailed re-analysis of the XMM-Newton EPIC observations of GRB 011211. In particular, 
we show that the detection of the soft X-ray line emission appears robust, regardless of detector background, 
calibration, spectral binning, or the spectral model that is assumed. We demonstrate that thermal emission, from 
an optically thin plasma, is the most plausible model that can account for the soft X-ray emission, which appears 
to be the case for at least two burst afterglow spectra observed by XMM-Newton. The X-ray spectrum of GRB 
011211 evolves with time over the first 12 ksec of the XMM-Newton observation, the observations suggest that 
thermal emission dominates the early afterglow spectrum, whilst a power-law component dominates the latter 
stages. Finally we estimate the mass of the ejected material in GRB 011211 to be of the order 4-20 solar masses. 
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1. Introduction 

The study of X-ray afterglow emission is crucial to under- 
standing the nature of gamma ray bursts and their pro- 
genitors, as X-ray spectroscopy can reveal details of the 
environment of the burst explosion. For instance, observa- 
tions of some X-ray afterglows with Beppo-SAX, ASCA, 
and Chandra have revealed strong iron Ka emission lines 
(e.g. GRB 000214, Antonelli et al. 2000; GRB 991216, 
Piro et al. 2000; GRB 970828, Yoshida et al. 1999; GRB 
970508, Piro et al. 1998), with line equivalent widths of 
up to several kcV; typically the lines are detected at ~ 3er 
confidence. An iron K-shell absorption feature was also re- 
ported in the prompt Beppo-SAX spectrum of one burst, 
GRB 990705 (Amati et al. 2000, Lazzati et al. 2001). 
The line observations appear to support models where 
the burst explodes into an enriched, high density medium, 
favoring a massive stellar progenitor for long duration 
bursts. In many cases (e.g. Piro et al. 2000, Antonelli et al. 
2000), large masses of iron are required to account for the 
line emission. One possible source of the iron is in a shell of 
distant (R ~ 10 16 cm) material ejected from a supernova 
(the 'supranova' model, Vietri & Stella 1998) which would 



have to occur at least several months prior to the burst 
event, in order to allow sufficient time for large quantities 
of iron to form (e.g. Vietri & Stella 1998). Alternatively, 
the iron K emission could arise through the interaction 
of a magnetically driven wind with the envelope of the 
massive progenitor star (Rees & Meszaros 2000); here the 
distances involved are much smaller (R^ 10 13 cm) and 
consequently the high masses of iron and long time delays 
between supernova and burst are not required. 

The high throughput of XMM-Newton, compared with 
Beppo-SAX or Chandra, makes it the best available tele- 
scope with which to constrain X-ray line emission in 
GRBs, potentially providing a powerful diagnostic for dis- 
criminating between models of the burst progenitor. In an 
analysis of the XMM-Newton observation of the afterglow 
of gamma ray burst GRB 011211, Reeves et al. (2002) 
reported the first detection of emission lines other than 
iron in a GRB afterglow spectrum; specifically the decay- 
ing line emission from the hydrogenic states of Mg, Si, S, 
Ar and Ca. Furthermore the energies of the X-ray lines 
appeared to be offset from the known redshift for the host 
galaxy of GRB 011211 (z = 2.140 ± 0.001, Holland et al. 
2002), implying the line emitting matter was outflowing 
with a velocity of ~ 0.1c from the site of the GRB pro- 
genitor. This result was interpreted as possible evidence of 
a supernova explosion occurring within days of the burst 
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itself, with the X-ray line emission arising from matter in 
the expanding supernova shell. 

Recently, Borozdin & Trudolyubov (2002) have 
claimed that the detection of soft X-ray emission lines 
from light metals in GRB 011211 may not be robust, 
as the observations could be contaminated by the back- 
ground of the XMM-Newton EPIC-pn detector. In addi- 
tion, Rutledge & Sako (2002) perform simulations which 
suggest that the significance of the emission features may 
be lower, at the level (at best) of ~ 98.8% confidence. In 
this paper, we describe our analysis of the XMM-Newton 
data on GRB 011211 in detail. We confirm that the soft X- 
ray emission features in GRB 011211 are robust, the set 
of lines being detected to a good level of confidence (at 
> 99% significance in all our tests). Both the F-test and 
Monte-Carlo simulations are used to determine the sta- 
tistical significance of the line emission features, making 
no prior assumption about the rest energies or redshift 
of the line emitting material. We also demonstrate that 
the results do not depend on the background subtraction 
of the EPIC-pn detector, instrument calibration, or the 
spectral binning used and show that the result is indepen- 
dent of the spectral model that is assumed. In section 2 
we review the properties of the burst, the XMM-Newton 
observations and basic data analysis and in section 3 de- 
tail the spectral analysis of the afterglow during the first 
5 ksec when the lines appear most prominent. The time 
dependent spectral properties of GRB 011211 are also out- 
lined, whilst a likely mechanism of the line emission is dis- 
cussed and a mass for the ejected material in GRB 011211 
is derived. Throughout this paper we adopt a cosmology 
of Ho = 75 km s _1 Mpc -1 and go = 0.1. Unless other- 
wise stated, all errors are quoted at 68% confidence (i.e. 
Ax 2 = 1.0 or 2.3 for 1 or 2 interesting parameters respec- 
tively) . 

2. XMM-Newton Observations of GRB 011211 

The Gamma Ray Burst GRB 011211 was first detected 
on December 11th 2001 at 19:09:21 (UT), by Beppo- 
SAX (Frontera et al. 2001). The burst duration was 
270s, making GRB 011211 the longest burst yet ob- 
served by Beppo-SAX, with a peak flux (40-700 keV) of 
5 x 10 -8 erg cm -2 s _1 . Spectroscopy of the transient opti- 
cal afterglow revealed several absorption lines at a redshift 
of z = 2.140±0.001 (Holland et al. 2002). Optical imaging 
(Burud et al. 2001, Fox 2002) has also linked the optical 
transient to extended emission, the probable host galaxy 
of magnitude m v = 25.0 ± 0.5. Assuming the absorption 
system arises from the host galaxy of the GRB, then the 
total equivalent isotropic luminosity for GRB 011211 is 
5 x 10 52 erg. 

The XMM-Newton (Jansen et al. 2001) observations 
of GRB 011211 started at 06:16:56 (UT) on December 
12th 2001, ~ 11 hours after the initial burst (Santos-Lleo 
et al. 2001). Data from EPIC (European Photon Imaging 
Cameras) have been analysed, using both the MOS CCDs 
(Turner et al. 2001) with the medium filter and pn CCDs 
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Fig. 1. EPIC 0.2-10keV lightcurve (coadded pn and 
MOS). A power-law decay fits the data well, the best- 
fit power-law, (with a decay index of 1.6 ± 0.1) is plotted. 
Time since the GRB is plotted on the abscissa. 

(Strudcr et al. 2002) using the thin filter. The total obser- 
vation duration (for EPIC-pn) was 27 ksec. Approximately 
4.2 ksec into the observation, the spacecraft pointing was 
moved from its original position at co-ordinates (2000) llh 
15' 16.4", -21°55' 44.8", to new boresight co-ordinates at 
llh 15' 18.0", -21°56' 53.0", in order to place the source 
away from the EPIC-pn chip gaps. The time-averaged 0.2- 
10 keV flux was 1.7 x 10~ 13 erg cm~ 2 s _1 (correspond- 
ing to an 0.6-30 keV, rest frame afterglow luminosity 
of 7 x 10 45 erg s _1 ), decreasing during the observation 
according to a power-law, with decay rate F(t) cx t~ a , 
a = 1.6 ± 0.1 (Fig. |). The X-ray decay is notably 
steeper than reported for the same period in the optical 
(a = 0.83 ± 0.04, Holland et al. 2002). 

2.1. Extraction of the XMM-Newton data 

We initially concentrate our analysis on the first 5 ksec of 
observation, where the prominent soft X-ray emission lines 
are observed. The analysis focuses on the EPIC-pn spec- 
trum, which due to the use of the thin filter, has a much 
higher count rate (> 6x) than the EPIC-MOS, having 
first checked that the MOS and pn data are self-consistent. 
All data reduction was performed using the XMM-SAS 
(v5.3). 

Fig. ||(a) shows the unscreened EPIC-pn field around 
GRB 011211, the afterglow itself being the bright source 
located close to the edge of the EPIC-pn chip gap, but near 
to the centre of the image. The background count rate is 
clearly higher around the edges of the pn chip in the un- 
screened dataset, as noted by Borozdin & Trudolyubov 
(2002), who claim that inclusion of the noisy pixels in the 
background spectrum could effect the detection of soft 
X-ray lines. However in order to produce a clean image 
and event list for GRB 011211, the dataset was screened 
to include only good X-ray events (using the selection 
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Fig. 2. XMM-Newton images of the field around GRB 
011211, during the first 5 ksec of observation; (a) top 
shows the unscreened image, whilst (b) bottom shows the 
screened image, used for extraction of all spectra. Here the 
bold circle indicates the source spectral extraction region 
used, of radius of 40". The dashed areas indicate the back- 
ground regions used in the extractions; (i) a box along the 
edge of the chip (but excluding the source region) or (ii) 
an offset circle of radius 80". 

expression 'FLAG=0', in the SAS tool evselect) using 
only the well-calibrated single or double pixel events (i.e. 
PATTERN<=4). An energy cut was also applied, exclud- 
ing all data below 200 eV and above 10 keV. The resulting 
image is shown in Fig. ^(b), which shows a smooth dis- 
tribution of background counts, regardless of position on 
the EPIC-pn chips. 

Source spectra for GRB 011211 were subsequently ex- 
tracted from a circular region of 40" radius, centred on 
the source, accounting for the decrease in flux incurred 
through photons falling onto the chip gaps. A background 
spectrum was taken from a circular region of radius 80", 
using a source-free region located close to the source 
(Fig. ||b). The background spectrum obtained is shown 
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Fig. 3. Background spectra used for spectral analysis dur- 
ing the initial 5 ksec, the grey circles show the spec- 
trum extracted from the offset circular background re- 
gion, whilst the crosses show the background spectrum ex- 
tracted from the box region, located close to the chip edge. 
The two background spectra are identical, and demon- 
strate that the background spectrum obtained is both fea- 
tureless and independent of the position on the pn chip. 
Overall the background count rate is ~ 10% of the net 
count rate from GRB 011211. 



in Fig. |3|; after applying the above screening the back- 
ground spectrum is smooth, unlike the background spec- 
trum obtained by Borozdin & Trudolyubov (2002), who 
using unscreened data obtain a 'spike' in their background 
spectrum near 0.7 keV. However, in order to ascertain 
whether the choice of background region can affect the 
spectrum of GRB 011211, we also use a background ex- 
tracted from a box-shaped area located close to the EPIC- 
pn chip gap (Fig. ^|b), but excluding the source region. The 
background spectra obtained for both regions are identi- 
cal (Fig. ||), hence the choice of background region does 
not effect the background subtracted source spectrum of 
GRB 011211. When normalized to the size of the source 
extraction region, the background contributes ~ 10% of 
the net source count rate. 

3. Spectral Analysis 

The background subtracted source spectrum of GRB 
011211, from the first 5 ksec of the observation, was fitted 
in XSPEC vll.l, using a minimum of 12 source counts per 
bin. An absorbing column of 4.2 x 10 20 cm~ 2 , due to our 
own galaxy, was included in all the spectral fits. A simple 
power-law fit (model 1) gives a barely (10% null hypoth- 
esis probability) acceptable fit to the data (\ 2 /d.o.f = 
54.2/45, where d.o.f is the number of degree of freedom), 
with a power-law photon index of T = 2.18 ± 0.10. The 
model fits to GRB 011211 are summarized in Table 1. 
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Fig. 4. The background subtracted XMM-Newton EPIC- 
pn spectrum of GRB 011211, from the first 5 ksec of ob- 
servation. Energy (units keV), in the observed frame, is 
plotted on the abscissa. In the upper panel, the EPIC-pn 
data points are plotted with filled circles, whilst the solid 
line shows a power-law model, folded through the EPIC- 
pn response function. The lower panel shows the ratio of 
the data to the power-law fit; an excess of counts is ob- 
served between 0.5 and 1.5 keV in the spectrum. This 
excess can be well fitted with at least 3 emission lines, 
at 0.71 ± 0.01 keV, 0.88 ± 0.01 keV and 1.21 ± 0.02 keV 
respectively. 



3.1. Soft X-ray Emission Line Fits to GRB 011211 

Inspection of the power-law fit to the spectrum of GRB 
011211 shows a clear excess present in the data residuals 
between 0.5 and 1.5 keV (Fig. ||). The deviations are con- 
siderably larger than the current calibration uncertainty 
of the EPIC-pn detector, which is of the order ~ 5%; 
adding a further systematic error of 5% to the data has 
no noticeable effect on the spectral fits. Hence, follow- 
ing the procedure of Reeves et al. (2002), we attempt 
to fit this excess with a set of emission lines. Initially 
3 narrow spectral lines were added (model 2), allowing 
the energy and flux of each line to vary, whilst the line 
widths were kept negligibly small (at a = 10 eV) in the 
fits, as the lines are unresolved compared to the detec- 
tor resolution. The best fit (observed frame) energies of 
the lines were then 0.71 ± 0.02 keV, 0.88 ± 0.01 keV and 
1.21 ± 0.02 keV, with line fluxes of (9.9 ± 3.4) x 10~ 6 pho- 
tons cm~ 2 s _1 , (7.7 ± 2.3) x 10~ 6 photons cm~ 2 s -1 and 
(4.0 ±1.4) x 10~ 6 photons cm -2 s -1 respectively. The fit is 
significantly improved, at 99.7% confidence, upon the ad- 
dition of these lines (x 2 /d.o.f = 33.7/39, A X 2 = 20.5 for 
6 extra degrees of freedom), when an F-test is employed 
(Bevington 1992). Note that even when fitted individu- 
ally, the lines are detected at 98%, 97% and 93% confi- 
dence (for Si, S and Ar respectively), with the addition of 
2 extra degrees of freedom (energy, flux) for each line. 

Examination of the data/model residuals in Fig. |J sug- 
gests that two further lines can be added to the spectrum 
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Fig. 5. The XMM-Newton spectrum of GRB 011211, fit- 
ted with the redshifted emission line model. The EPIC-pn 
spectrum is shown in black whilst the lower signal to noise 
EPIC-MOS spectrum is shown in grey. The lines fitted to 
the data correspond to Mgxi (or Mgxn), Sixiv, Sxvi, 
Ar XVIII and Caxx at a redshift of z = 1.90 ±0.06, imply- 
ing an outflow velocity of ~ 0.1c for the emitting material. 

Table 1. Summary of X-ray spectral fits to the first 5 
ksec of observation for GRB 011211. a Fit statistic, in 
terms of x 2 divided by degrees of freedom. b Significance 
of fit improvement, compared to a power-law, from using 
an F-tcst. ' indicates parameter is fixed in fit. 



Model 


T or kT 


X Z /d.o.f. a 


Prob" 


1. 


PL Only 


2.18 ±0.10 


54.2/45 




2. 


PL + 3 lines 


2.15 ± 0.11 


33.7/39 


99.7% 


3. 


PL + 5 lines 


2.10 ±0.15 


29.0/35 


99.3% 


4. 


PL + z-lines 


2.12 ± 0.12 


34.3/39 


99.5% 


5. 


PL + VMEKAL 


r = 2 s 


38.1/41 


99.8% 






kT = 4.1 ± 1.6 







(model 3), with observed energies 0.44 ±0.04 keV, 1.46 ± 
0.04 keV and fluxes (1.0 ± 0.8) x 10~ 5 photons cm~ 2 s" 1 
and (1.9 ± 1.1) x 10~ 6 photons cm~ 2 s^ 1 respectively, al- 
though the improvement in fit statistic is marginal for 
these last two cases (x 2 /d.o.f = 29.0/35). Table 2 lists 
the properties of each of the individual lines added to the 
X-ray spectra. Fig. | shows the EPIC-pn and EPIC-MOS 
data fitted with the 5 emission lines. Unfortunately none 
of the lines are significantly detected in the MOS CCDs, 
due to the lower count rate, although the data for both 
detectors are consistent (e.g. see the residuals in the lower 
panel) . 

When converted to the rest frame of the burst (at 
z=2.14), the energies of all 5 possible lines are 1.38 ± 
0.10 keV, 2.22±0.05keV, 2.77±0.05 keV, 3.80±0.06keV 
and 4.58 ± 0.12 keV. Interestingly these lines all appear 
to be offset from the known energies of the strong K-shell 
lines of Mgxi (1.35 keV), Sixiv (1.99 keV), Sxvi (2.60 
keV), Ar xvni (3.30 keV) and Caxx (4.07 keV). With the 
exception of Mgxi (the most tentative detection), these 
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Table 2. Soft X-ray line fits to GRB 011211. ° Known (laboratory) energy of emission line (in keV); b observed energy 
(in keV) of emission line in spectrum; c measured rest-frame (at z=2.14) energy (keV) of the emission line; d redshift 
for the emission line, calculated from a and b . Observed line flux, in units 10 -6 photons cm~ 2 s _1 . 



Ident 


77 a 
Inline 


Eobs 6 


E t c 

^rest 


redshift" 


Flux e 


Mgxi 


1.35 


0.44 ± 0.04 


1.38 ±0.10 


2.03 ± 0.22 


10 ±8 


Si xiv 


1.99 


0.71 ±0.02 


2.22 ±0.05 


1.82 ± 0.06 


9.9 ±3.4 


Sxvi 


2.60 


0.88 ±0.01 


2.77 ±0.05 


1.94 ±0.02 


7.7 ±2.3 


Ar xvni 


3.30 


1.21 ±0.02 


3.80 ±0.06 


1.73 ±0.05 


4.0 ± 1.4 


Caxx 


4.07 


1.46 ±0.04 


4.58 ±0.12 


1.79 ±0.07 


1.9 ± 1.1 



lines appear to correspond to the Lyman-a transitions of 
the most abundant elements, and are the strongest emis- 
sion lines one would expect to observe at these X-ray en- 
ergies for either a collisionally ionised plasma (of temper- 
ature a few keV) or a highly photoionised plasma. 

Having established the likely identity of the soft X-ray 
line features in GRB 011211, we then calculated the ob- 
served redshift, and thus the outflow velocity of the line 
emitting material. Fitting the redshift of the lines indi- 
vidually, one finds values of z = 2.03 ± 0.22 (Mgxi), 
z = 1.82±0.06 (Sixiv), z = 1.94±0.05 (Sxvi), 1.73±0.05 
(Arxvni) and 1.79 ± 0.07 (Caxx). Thus with the excep- 
tion of Mgxi (which at the detector resolution is unre- 
solved from Mgxn and also L-shell emission from Fe or 
Ni) the lines all appear to be significantly offset from the 
known redshift for GRB 011211 of z = 2.140 ± 0.001, cor- 
responding to an outflow velocity of ~ 0.1c. In order to 
test the significance of this, we fixed the lines at these 
known rest-frame energies, but allowed a single redshift 
(between z — and z = 10, in steps Az = 0.05) for the 
entire line set to vary (model 4). The redshift calculated 
for the line set is z = 1.90 ± 0.06, corresponding to an 
outflow velocity of 0.085 ± 0.02c or - 25000 km s" 1 from 
the GRB progenitor. Overall, the fit statistic improves to 
X 2 /d.o.f = 34.3/39. When compared to the power-law fit 
(for 6 extra degrees of freedom), the addition of the red- 
shifted line set is significant at 99.5% confidence from an 
F-test. 



3.2. Thermal Emission and Ionised Reflection models 

Finally, we attempted to fit the data with emission from a 
collisionally ionised optically thin plasma, with variable 
elemental abundances (hereafter the 'thermal' emission 
model), using the vmekal code within xspec (Mewe et 
al. 1985). Initially, a thermal emission component of tem- 
perature kT = 4.1 ± 1.6 keV was added in addition to an 
underlying power-law component (model 5) . The improve- 
ment in fit obtained (A^ 2 /d.o.f — 38.1/41) is significant 
when compared to the pure power-law fit (model 1); the 
thermal emission component is required at 99.8% confi- 
dence from an F-test, for four extra degrees of freedom 
(temperature, abundance, redshift and luminosity). Note 
that the spectrum during the initial 5 ksec can be equally 
well fitted with pure thermal emission (i.e. the power-law 
component is not statistically required) the fit statistic is 
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Ar XVIII 
Ar XVII Ca XX 

Ca XIX 



2 3 4 

Rest Energy (keV) 

Fig. 6. Rest-frame thermal emission model showing the 
relative strengths of the major emission lines (and contin- 
uum) for a plasma of temperature 4 keV. 

then x 2 /d.o.f = 38.5/43. The rest-frame model spectrum 
using the vmekal component is shown in Fig. ^. It can be 
seen that for these elements, the hydrogenic states domi- 
nate (e.g. Si XIV at 2.0 keV), although there is also some 
contribution from the He-like species as well (e.g. Sxv, 
Arxvn, Caxix). 

We also used this model to provide an independent de- 
termination of the redshift of the line emitting material, 
where no assumption has been made about the likely iden- 
tity of the lines; the best-fit redshift was 1.89 ± 0.08, in 
excellent agreement with the simple line fit derived above. 
A contour plot showing the 68%, 90%, 95% and 99% con- 
fidence levels for the line emitting material is shown in 
Fig. 0, a redshift of z = 2.14 can be ruled out at > 99% 
confidence. Rutledge & Sako (2002) suggest that emission 
at z=1.2 and z=2.75 are also plausible solutions to the 
data, clearly these are ruled out from our fits, although 
there is a local minimum at z ~ 1.2 (which can be ex- 
cluded at ~ 99% confidence). Nonetheless, we consider 
both the lower redshift of z = 1.2 and the high redshift 
value of z = 2.75 to be physically implausible, as this 
would require either an extreme blueshift or for the GRB 
to reside at a higher redshift than the measured value 
(from the optical spectrum) of z = 2.14. 

The relative abundance of the light metals (Mg, Si, 
S, Ar, Ca) was found to be ~ 9x the solar value, with 
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Fig. 7. Confidence contour plot of measured flux against 
redshift for the thermal line emission model, added in ad- 
dition to a power-law component. The contours represent 
the 68%, 90%, 95%, 99% confidence levels to the spectral 
fit, whilst the redshift has been varied (in Az = 0.05 incre- 
ments) between z = 1 and z — 3. The thermal emission 
is required in the data at > 99% confidence, whilst the 
need for an outflow velocity, with respect to the GRB rest 
frame, is also required at > 99% confidence. 



a 90% lower-limit of 4x solar. Interestingly, unlike the 
other elements, an over abundance of iron is not required, 
the 90% upper limit for the iron abundance was found 
to be < 1.4x solar, consistent with the lack of detection 
(EW< 400 eV, rest frame) of a strong iron K-shell line in 
the data. 

We also attempted to fit the spectrum with a reflec- 
tion model (Ballantyne & Ramirez- Ruiz 2001), where the 
X-rays scatter off the walls of a funnel excavated in the 
envelope of a stellar progenitor (the so-called "nearby 
reprocessor" model, Rees & Meszaros 2000). Firstly, we 
tried a pure reflection component (i.e. just the scattered 
X-ray emission). Allowing the ionisation state of the re- 
flecting material to vary results in an ionisation param- 
eter log £ > 4 (with r = 2.0) and a fit statistic of 
X 2 /d.o.f — 55.3/44. In this fit, the spectrum resembles 
a simple power-law, as all the elements are fully ionised, 
and the reflector is a perfect mirror. If one forces the ion- 
isation parameter to be log £ < 4, then the fit statistic 
is worse (x 2 /d.o.f — 66.7/43), whilst the incident power- 
law becomes steep (r ~ 3). The poor-fit results from the 
strong iron line or edge present in the model, and from the 
shape of the continuum, which is concave in the reflection 
model (but is convex in the actual data). Allowing a di- 
rect power-law component in addition to the reflector does 
not improve the fit further. Hence it is difficult to reconcile 
the current reflection models with the spectrum of GRB 
011211, although we note that at present, this model does 
not include emission from S, Ar or Ca. 



3.3. Monte-Carlo Simulations 

We have used Monte-Carlo techniques to provide an inde- 
pendent estimate of the significance of the line emission 
fits to GRB 011211. This is important given the low pho- 
ton number statistics of the data, where the distribution of 
counts may become non-Gaussian (Protassov et al. 2002). 

In order to test the significance of the line emission, 
10000 random spectra from the EPIC-pn detector were 
generated, using a simple power-law emission model (in- 
cluding absorption from our own Galaxy), with a photon 
index of T — 2.2. An exposure time of 5 ksec (accounting 
for detector deadtime) and an X-ray flux identical to that 
of GRB 011211 (F . 2 _i 0fce y = 2.1 x 10~ 13 erg cm" 2 s" 1 ) 
were used, as well as a background level ~ 10% of the 
source count rate, so that the signal to noise of the simu- 
lated spectra were comparable to the actual source spec- 
trum for GRB 011211. The simulated spectra were re- 
binned to obtain a minimum of 12 source counts in each 
spectral channel. 

We proceeded to fit the simulated data, using (i) a 
pure power-law continuum and (ii) a model consisting 
of a power-law plus three emission lines, corresponding 
to the strongest lines (Si, S, Ar) observed in the actual 
GRB 011211 spectrum. The difference in Ax 2 between 
the power-law and line emission models was noted for each 
simulated spectrum, in order to test whether the appar- 
ent line emission could be produced purely through ran- 
dom, Poisson noise. Intially, the line energies were fixed 
at the rest-frame values corresponding to the Ka tran- 
sistions of Sixiv, Sxvi and Arxvni, but note that the 
redshift was allowed to vary between z = and z = 10 
(stepsize Az = 0.05) when fitting the trial spectra. Hence 
no prior assumption is made about the redshift of the line 
emitting material. In addition, the line flux for each line 
was allowed to vary in the fits. Upon fitting the simulations 
we found that only 5 in 10000 of the randomly generated 
power-law spectra could yield a similar decrease in chi- 
squared (A% 2 ~ 15), when compared to the real data, af- 
ter the addition of the three emission lines. The likelihood 
(or null hypothesis probability) that the lines could arise 
purely from Poisson noise in this test is < 0.05%; hence 
the simulations indicate that the soft X-ray line emission 
is detected at 99.95% confidence, in agreement with the 
results obtained through the F-test. 

In a recent paper, Rutledge & Sako (2002) note that 
the prior identification of the three strongest emission line 
features with Sixiv, Sxvi and ArXVlll represents an a 
posteriori assumption about the data, when trying to esti- 
mate the significance of the line features. However we note 
(and see Fig. ||) that line emission features from these 3 
abundant H-like species will dominate the emission from a 
hot (several keV) collisionally ionised plasma, or a highly 
photoionised plasma, over the 1-4 keV (rest frame) en- 
ergy range. Nonetheless we relax the assumption that the 
lines originate from these species, and test the improve- 
ment in fit obtained when adding lines, at any energy from 
0.2-10 keV in the observed spectrum, to the 10000 fake 
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Fig. 8. Histograms showing the improvements in the 
fit statistic (x 2 ) obtained when fitting 10000 simulated 
power-law spectra with; (a) a model consisting of a power- 
law plus three emission lines (with the line energies 
and fluxes both allowed to vary) and (b) the thermal 
(vmekal) emission model. In both cases, only three out 
of 10000 randomly generated trial spectra yield an im- 
provement in x 2 equivalent to that found in the actual 
observation; note the A% 2 for the real GRB 011211 data is 
marked by an arrow. The simulations show that the prob- 
ability of generating the soft X-ray line emission through 
random Poisson noise is only 0.03%. 



power-law spectra (i.e. a blind spectral search). Initially 
we just add one line to the fake spectra (whilst varying 
both the energy and flux of the line), and find that an im- 
provement of A% 2 = 8, equivalent to the addition of the 
Si xiv line in the real dataset, is obtained in only 6 out of 
1000 cases. Thus the Sixiv line is detected at 99.4% con- 
fidence. Similarly we find that the S xvi and Ar xvm lines 
are detected at 98.8% and 96% confidence respectively, in 
agreement with the values computed by the F-test. We 
also consider the significance of the 3 line set as a whole, 
by fitting 3 Gaussian lines (again at arbitary energies) to 
the simulated datasets. The results (see Fig. ||a) show that 
only 3 in 10000 trial spectra gave an equivalent (or better) 
improvement in the \ 2 statistic compared to the real data 
(where A% 2 = 20.5 for the addition of all 3 lines). Thus 
we can conclude that the addition of 3 Gaussian lines in 
the soft X-ray spectrum of GRB 011211 is significant at 
> 99.9% confidence, whilst the individual lines are de- 
tected at 99.4%, 98.6% and 96% confidence (for Si, S and 
Ar) respectively. 

A more physically realistic model, consisting of emis- 
sion from a hot, thermal plasma (i.e. the vmekal model), 
was also fitted to the 10000 randomly generated power-law 
spectra. The temperature, luminosity, redshift and rela- 
tive abundance of the plasma were allowed to vary in the 
fits. Fig. H (panel b) shows a histogram of the change in 
\ 2 obtained when comparing a power-law fit to a thermal 
fit for the simulated spectra; note here a positive change 



in x 2 corresponds to an improvement in fit statistic when 
comparing the vmekal model to a power-law. In the ma- 
jority (90%) of cases, fitting a thermal model to the ran- 
domly generated power-law spectra results in a worse fit. 
However the real data for GRB 011211 indicate that an 
improvement in fit statistic of A% 2 = 15.7 is obtained us- 
ing the thermal model compared to the pure power-law. 
Indeed an equivalent improvement in A% 2 is only found 
in 3 out of 10000 simulated spectra; hence the inclusion of 
the thermal line emission component is formally required 
at > 99.9% confidence from the simulations. 

3.4. Other Extraction Effects 

We have also investigated whether varying the size of the 
source selection region has any effect on the derived spec- 
trum and the detection of the soft X-ray lines. In order 
to test this, we extracted spectra from circular regions of 
radii 10", 20" 30" in addition to the previous spectrum 
that was extracted with a 40" radius. We find the detec- 
tion of the redshifted soft X-ray line set is significant in 
each case; at >99% confidence (via an F-test) for spectra 
extracted with 20" and 30" radii, whilst the detection is 
less significant (at 98% confidence) for a 10" radius (as 
substantially fewer source counts are included in the spec- 
trum). However, after accounting for the telescope point 
spread function, the derived line flux for the Si Xiv line was 
found to be constant in each case (at ~ 1.0 x 10~ 5 pho- 
tons cm~ 2 s _1 ), in contrast to the findings of Borozdin & 
Trudolyubov (2002) who claim the Si xiv line flux in their 
spectra decreases with the extraction radius. 

The effect of spectral binning was also investigated. 
Rather than binning the spectrum to a set number of pho- 
tons per bin, the spectra were re-binned so that each bin 
represented a constant energy interval, of 40 eV in size. 
This is then well matched to the resolution of EPIC-pn 
detector below 2 keV (where the FWHM spectral resolu- 
tion is ~ 80 eV). Data above 2 keV were ignored due to 
poor photon statistics (less than 10 source counts per bin). 
The resultant spectrum is shown in Fig. ^|. The spectrum 
has been fitted with the 5 redshifted lines obtained pre- 
viously, the fit statistic again improves significantly (from 
x 2 /dof = 60.5/42 to 33.4/32, i.e. at 99% confidence) upon 
the addition of the line set. Interestingly, a residual near 
1.1 keV is present in the data, the addition of a sixth 
emission line, at 1.09 ± 0.02 keV, is significant at 98% 
confidence (x 2 /d.o.f = 25.9/30). Assuming a redshift of 
z = 1.9, this line could correspond to He- like Ar (ArXVll), 
which is predicted in the thermal model (Fig. |^). The net 
redshift obtained for the line set is then z — 1.86 ± 0.04, 
consistent with the previous determination (section 3.1). 

4. Temporal variability of the Line Emission in 
GRB 011211 

The observations presented in Reeves et al. (2002) showed 
that the line emission in GRB 011211 decayed with time 
faster than the continuum, rendering the lines unde- 
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Fig. 9. EPIC-pn spectrum of GRB 011211, binned into 
equal energy bins. 40 eV in size. The 5 redshift emission 
lines are clearly detected, and a sixth line at 1.1 keV, due 
to He-like Ar (Arxvn) may also be present (at 98% con- 
fidence) . The likely identifications of the lines are marked 
at the top of the plot. 



tectable after 10 ks of observation. To investigate this fur- 
ther, we initially split the first 12 ksec of observations into 
4x3 ksec segments. Although it is difficult to constrain the 
significance of the individual lines over such a short time 
period, nonetheless we compare the power-law model with 
the thermal (vmekal) model and the model consisting of 
a power-law plus rcdshiftcd emission lines. The results of 
the fits are shown in table 3. 

It can be seen that the two spectra from the initial 
6 ksec of observation are poorly fit by a simple power- 
law spectrum, whilst both the redshifted line and thermal 
emission models are a significantly better fit. What ap- 
pears interesting is that the early afterglow spectrum is 
significantly curved. The spectrum during the first 3 ksec 
can be fit with a steep power-law (with T = 3.6 ± 0.7) 
and an observed frame column density of Nh = 3.0 ± 
1.1 x 10 21 cm~ 2 . This indicates that the early afterglow 
spectrum is dominated by an intrinsically curved (ther- 
mal) component. In contrast, acceptable fits to the last 
two spectra (between 6 and 12 ksec) can be obtained by 
a pure power-law model, neither the soft X-ray line emis- 
sion, nor the addition of a thermal emission component 
being formally required in these spectra. 

In order to quantify this, we have calculated the ra- 
tio of thermal to power-law emission during the first 
and last 10 ks of observation. During the first 10 ksec, 
the thermal component has an average luminosity of 
2.0 x 10 45 erg s _1 , whilst the power-law component has a 
luminosity of 6.8 x 10 45 erg s -1 (and note that both compo- 
nents are statistically required in the spectral fit at > 99% 
confidence). In the last 10 ksec the power- law component 
has a luminosity 4.85 x 10 45 erg s _1 , whilst the thermal 
emission is no longer required (L < 1.07 x 10 45 erg s" 1 ). 



We have also parameterised the strength of the line 
emission by measuring the equivalent width (EW) of the 
Si xiv line, relative to a pure power-law continuum. During 
the first 5 ksec, the Sixiv line is detected with an (ob- 
served frame) EW of 127±42 eV. A line is also marginally 
detected (at 90% confidence) from 5-10 ksec after the start 
of the observation, with a smaller EW of 62 ± 30 eV. 
Beyond 10 ks, the line is not detected, we can place a 
tight upper-limit of < 25 eV on the EW of the line from 
the last 17 ks of observation. Thus clearly the soft X-ray 
line emission in GRB 011211 decays with time, and is not 
detected after 10 ks from the start of the XMM-Newton 
observation. 

5. Discussion and Conclusions 

Both the F-test and Monte-Carlo simulations have showed 
that the soft X-ray line emission in GRB 011211 is de- 
tected at > 99% confidence, We conclude that the detec- 
tion of the soft X-ray features in GRB 01211 is robust, 
and is not affected by the detector background, calibra- 
tion, spectral binning or the particular spectral model that 
is assumed. It is also apparent that the lines decay rapidly, 
no line emission is formally detected after 10 ks from the 
start of the XMM-Newton observation. 

The early afterglow data (during the first 10 ksec of 
observation) appears to be characterised by thermal X- 
ray emission (from optically thin gas), as evident by the 
curvature present in the in the early spectrum. An alter- 
native spectral model, where the afterglow line emission 
results from X-ray reflection off optically thick material 
appears to be ruled out (e.g. Ballantyne & Ramirez-Ruiz 
2001), although we note that in current reflection models, 
the emission from S, Ar and Ca is not computed. Future 
models will hopefully provide a better test of the reflection 
scenario. 

On the assumption that the lines result from thermal 
emission one can estimate the mass of the ejected material. 
Firstly the emission measure (Em) is computed from the 
observed luminosity of the thermal emission component. 
Here n 2 V = Em ~ 10 69 cm~ 3 (where n is the electron den- 
sity in cm -3 and V is the volume), whilst the plasma will 
emit X-rays over a characteristic cooling time (governed 
by the electron density) of t coo i = 1.4 x 10 15 n _1 sees, for 
gas of temperature kT = 4 keV. Thus from the emission 
measure, the total integrated output of the thermal emis- 
sion is governed by n 2 V = E M tii ne /t coo u where t Une is the 
observed lifetime of the line emission. Substituting for the 
cooling time, and as nV = M/m p (where M is the total 
mass of the emitting material and m p is the proton mass) , 
then one can derive the total mass of the ejecta:- 

M = 2.4 x 10 30 t ime g 
or M = 1.2 x 10- 3 ^ me Mq 

The minimum lifetime of the line emission (mea- 
sured from the start of the XMM-Newton observation) 
is 10 4 /(1 + z) sees, whilst the maximum lifetime occurs 
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Table 3. Fits to GRB 011211, over the first 12 kscc of observation. a Time measured from start of observation. b Power- 
law model with Galactic absorption. c Power-law model with free neutral absorption (measured in observed frame). 
d Units of column density 10 21 cm" 2 . e Fit using variable abundance thermal plasma model vmekal. ' Indicates 
parameter is fixed in fit. 9 Rcdshift determined from model. h Power-law model including 5 rcdshiftcd emission lines 
from Mg, Si, S, Ar and Ca. 



Model 




0-3 ksec" 


3-6 ksec" 


6-9 ksec a 


9-12 kscc a 


PL only" 


r 


2.15 ±0.16 


2.32 ±0.21 


2.14 ±0.06 


2.42 ± 0.20 




X 2 /d.o.f 


51.5/28 


23.9/16 


26.8/24 


18.7/26 


PL*wabs c 


r 


3.6 ±0.7 


4.3 ± 1.1 


2.9 ±0.5 


2.42 ± 0.20 




N H d 


3.0 ± 1.1 


3.6 ± 1.5 


1.6 ±0.8 


< 2.5 




X 2 /d.o.f 


39.1/27 


12.1/15 


23.6/23 


16.9/26 


Thermal 6 


kT 


4.1 ± 1.2 


3.3 ±0.8 


3.1 ±0.6 


2.3 ±0.5 




■z? 


1.89 ±0.09 


1.89 ±0.12 


1.9 / 


1.9 / 




X 2 /d.o.f 


31.9/26 


14.0/14 


20.6/22 


14.1/24 


PL + z-lincs' 1 


r 


3.0 ± 1.1 


2.4 ±0.5 








z 9 


1.90 ±0.08 


1.86 ±0.07 








X 2 /d.o.f 


21.9/22 


11.4/10 







if the line emission starts at the time of the burst, which 
is then time = 5.5 x 10 4 /(1 ± z) sees. Thus the minimum 
and maximum mass of the outflowing material, responsi- 
ble for the blue-shifted line emission, are then 4 and 20 
solar masses respectively. This indicates that the mass of 
the stellar progenitor is likely to be > 20M Q even for the 
conservative case where the mass of the ejecta is 4M Q (e.g. 
Woosley & Weaver 1995). 

The duration of the line features in GRB 011211, ob- 
served 11 hours after the initial burst, implies a distance 
of several x 10 15 cm for the line emitting matter (Reeves 
et al. 2002). When combined with the velocity of the ma- 
terial (~ 0.1c), this implies that the matter was ejected 
from days to weeks before the burst itself occurred. This 
time delay appears too short in the context of the "supra- 
nova" model (Vietri & Stella 1998), where a time-delay of 
several months to years can be predicted, whilst the "hy- 
pernova" class of models (Woosley 1993) involve the near- 
simultaneous occurrence of a supernova and the gamma 
ray burst. Interestingly one model has been proposed in- 
volving the coalescence of a close binary system, resulting 
from a collapse of a massive stellar progenitor (Davies et 
al. 2002). This could naturally account for the time delay 
between the putative supernova and the gamma ray burst, 
with a delay of a few days resulting from the orbital de- 
cay (through gravitational radiation) of the close binary 
system. 

If the line emitting material is located at much smaller 
distances (R ~ 10 13 cm), as in nearby reprocessor models 
(e.g. Rees & Meszaros 2000) then the need for a signifi- 
cant time delay is reduced. Here, the line emission results 
from re-processing (via reflection/photoionisation) of con- 
tinuum photons emitted after the initial burst. However in 
this scenario, significant continuum emission should be ob- 
served at the same time as the line emission, which would 
render any soft X-ray line features undetectable by cur- 
rent instrumentation (see Lazzati et al. 2002). One possi- 
bility that has been suggested (Kumar & Narayan 2002) is 
that the ionising photons are scattered back towards the 



ejected material, by a positron-electron shield generated 
by the initial 7-rays themselves. This model removes the 
need for a supernova-GRB delay, as the duration of the 
line emission is accounted for by the time required for the 
photons to be scattered back towards the ejected matter. 

Including GRB 011211, there have now been 4 de- 
tections of GRB afterglows with XMM-Newton (GRB 
001025A and GRB 010220, Watson et al. 2002; GRB 
020322, Ehle et al. 2002), whilst the X-ray afterglow 
for the faint burst GRB 020321 could not be localised. 
Interestingly, line emission features are also indicated in 
the X-ray afterglows of GRB 00 1025 A and GRB 010220 
(Watson et al. 2002). Similar to GRB 011211, the spec- 
trum of GRB 00 1025 A appears to exhibit a blend of soft 
X-ray line emission, from the medium-z elements such as 
Mg, Si, S, and Ar. In contrast, the spectrum of GRB 
010220 shows a strong iron-group (i.e. Fe, Co or Ni) line 
(equivalent width ~ 1 keV), although here the X-ray spec- 
trum is absorbed below 1 keV, rendering any soft X-ray 
line emission undetectable. Both fits favour an over abun- 
dance of Nickel (or Cobalt) to iron, which is suggestive of 
a short time delay (days rather than months) between the 
putative supernova and burst. Thermal emission models 
are formally required to fit the soft X-ray excess observed 
in GRB 001025A, whilst the spectrum of GRB 010220 can 
be equally well fitted by either thermal or reflection mod- 
els. Thus soft X-ray line emission is required in at least 
2 of the XMM-Newton afterglow spectra (GRB 011211 
and GRB 00 1025 A), the exception being the recent burst 
GRB 020322, which appears to exhibit a featureless (but 
absorbed) power-law spectrum (Reeves et al. 2002b, in 
preparation) . 

It seems possible that soft X-ray line emission features 
are relatively common in X-ray afterglow spectra, whilst 
thermal emission models can fit 3 out of the 4 present 
XMM-Newton spectra. The previous lack of detection of 
soft X-ray features with Chandra and Beppo-SAX can be 
accounted for by the lower effective area of the ACIS-S 
and LECS detectors respectively in the soft X-ray band. 
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Indeed XMM-Newton is the first X-ray mission with the 
sensitivity required to detect any putative soft X-ray fea- 
tures in GRBs at energies below 2 keV. It is important that 
XMM-Newton performs follow-up observations of bright 
(> 1CU 13 erg cm~ 2 s _1 ) afterglows, within several hours 
of the burst, to determine the frequency the soft X-ray 
line emission. Ultimately, through monitoring the tempo- 
ral behavior of the line emission from minutes to many 
hours after a burst, the X-ray telescope on-board NASA's 
forthcoming Swift mission will be able to trace the distri- 
bution of matter around the burst explosion and solve the 
ambiguities that arise when discriminating between the 
current GRB models. 
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